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ABSTRACT 

The improved Cascade-Exciton Model code CEM2k 
and the Los Alamos version of the Quark-Gluon 
String Model code LAQGSM, previously merged with 
the Generalized Evaporation Model code of Furihata 
(GEM2) were further modified to provide reliable 
proton- induced fission cross sections for applications. 
By adjusting two parameters in GEM2 for each mea- 
sured reaction, we were able to describe very well with 
CEM2k+GEM2 and LAQGSM+GEM2 all available ex- 
perimental fission cross sections induced by protons 
with energies from 20 MeV to 10 GeV both for subac- 
tinide and actinide targets. We also successfully tested 
our approach on several reactions induced by neutrons, 
pions, and photons. 



Introduction 

In recent years, an improved version of the Cascade- 
Exciton Model (CEM), contained in the code CEM2k 
[1] and the Los Alamos version of the Quark-Gluon 
String Model, implemented in the high-energy code 
LAQGSM [2] have been developed at the Los Alamos 
National Laboratory for a number of applications. 
CEM2k is intended to describe nucleon-, pion-, and 
photon-induced reactions at incident energies up to 
about 5 GeV, while LAQGSM describes both particle- 
and nucleus-nucleus reactions at energies up to about 1 
TeV/nucleon. Originally, both CEM2k and LAQGSM 
were not able to describe fission reactions and pro- 
duction of light fragments heavier than 4 He, as they 
had neither a high-energy-fission nor a fragmentation 
model. Recently, we addressed these problems [3, 4] 
by further improving our codes and by merging them 
with the Generalized Evaporation Model code GEM2 
developed by Furihata [5]. 

GEM2 is an extension by Furihata of the Dostrovsky 
et al. [6] evaporation model as implemented in LAHET 



[7] to include up to 66 types of particles and fragments 
that can be evaporated from an excited compound nu- 
cleus plus a modification of the version of Atchison's fis- 
sion model [8] used in LAHET. It was found [3, 4] that 
if we were to merge GEM2 with CEM2k or LAQGSM 
without any modifications, the new code would not de- 
scribe correctly the fission cross section (and the yields 
of fission fragments). This is because Atchison fitted 
the parameters of his fission model when it was cou- 
pled with the Bertini fntra-Nuclear Cascade (INC) [9] 
which differs from our INC. In addition, Atchison did 
not model prccquilibrium emission. Therefore, the dis- 
tributions of fissioning nuclei in A, Z, and excitation 
energy E* simulated by Atchison differ significantly of 
the distributions we get; as a consequence, all the fission 
characteristics are also different. Furihata used GEM2 
coupled either with the Bertini INC [9] or with the IS- 
ABEL [10] INC code, which also differs from our INC, 
and did not include preequilibrium particle emission. 
Therefore the real fissioning nuclei simulated by Furi- 
hata differ from the ones in our simulations, and the 
parameters adjusted by Furihata to work the best with 
her INC will not be the best for us. To get a good de- 
scription of fission cross sections (and fission-fragment 
yields) we need to modify at least two parameters in 
GEM2 (see details in [3, 4]). This problem was solved 
both for CEM2k+GEM2 and LAQGSM+GEM2 in the 
present work. 



Calculation of a f in GEM2 

A comprehensive description of GEM2 was published 
by Furihata [5], some details may be found in our pa- 
pers [3, 4], therefore we recall here only how fission 
cross sections are calculated by GEM2, as we need to 
modify them here. The fission model used in GEM2 is 
based on Atchison's model [8] , often referred in the lit- 
erature as the Rutherford Appleton Laboratory (RAL) 
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model, which is where Atchison developed it. There are 
two choices of parameters for the fission model: one of 
them is the original parameter set by Atchison [8] as 
implemented in LAHET [7] , and the other is a parame- 
ter set evaluated by Furihata [5] , used here as a default 
of GEM2. 

The Atchison fission model is designed to only de- 
scribe fission of nuclei with Z > 70 (we extended it 
in our codes down to Z > 65). It assumes that fis- 
sion competes only with neutron emission, i.e., from 
the widths Tj of n, p, d, t, 3 He, and 4 He, the RAL code 
calculates the probability of evaporation of any particle. 
When a charged particle is selected to be evaporated, 
no fission competition is taken into account. When a 
neutron is selected to be evaporated, the code does not 
actually simulate its evaporation, instead it considers 
that fission may compete, and chooses either fission or 
evaporation of a neutron according to the fission prob- 
ability Pf. This quantity is treated by the RAL code 
differently for the elements above and below Z — 89. 

1) 70 < Zj < 88. For fissioning nuclei with 70 < 
Zj < 88, GEM2 uses the original Atchison calculation 
of the neutron emission width T n and fission width T f 
to estimate the fission probability as 



Pf 



i 



r / + r„ i + r„/iy 



(i) 



Atchison uses [8] the Weisskopf and Ewing statistical 
model [11] with an energy-independent pre-exponential 
factor for the level density and Dostrovsky's [6] inverse 
cross section for neutrons and estimates the neutron 
width T„ as 

T„ = 0.352(l.68J + 1.93^ 1/3 Ji 
+^ /3 (0.76Ji -0.05 Jo)), 

where Jo and J\ are functions of the level density pa- 
rameter a n and s n (— 2y 'a n (E — Q n — 5)) as 



Jo 



(s n - l)e s " + 1 



Ji = 



2a n 

(2s 2 n - 6s n + 6)e s " + 4 - 6 
Sal 



The RAL model uses a fixed value for the level density 
parameter a n , namely 

a n = {Ai - l)/8. 

The fission width for nuclei with 70 < Zj < 88 is cal- 
culated in the RAL model and in GEM2 as 

(Sf - l)e°f + 1 

1 / = , 

a/ 



where Sf — 2^Jaf{E — Bf — 8) and the level density 
parameter in the fission mode a/ is fitted by Atchison 
to describe the measured V f/T n as 



a f = a„(l.08926 + 0.01098(x - 31.08551) 2 ) , 
and X = Z 2 /A. 



(2) 



2) Zj > 89. 



For heavy fissioning nuclei with Zj > 89 



GEM2 follows the RAL model and does not calculate 
at all the fission width Tf and does not use Eq. (1) to 
estimate the fission probability Pf. Instead, the follow- 
ing semi-empirical expression obtained by Atchison by 
approximating the experimental values of r„/r/ pub- 
lished by Vandenbosch and Huizenga [12] is used to 
calculate the fission probability: 



log(T n /T f ) = C(Z i )(A i -A (Z i )), 



(3) 



where C(Z) and Aq(Z) are constants dependent on the 
nuclear charge Z only. The values of these constants 
are those used in the current version of LAHET [7] and 
are tabulated in Table 1 (note that some adjustments 
of these values have been done since Atchison's papers 
[8] were published). 

Table 1. C(Z) and A (Z) values used in GEM2 



z 


C{Z) 


A (Z) 


89 


0.23000 


219.40 


90 


0.23300 


226.90 


91 


0.12225 


229.75 


92 


0.14727 


234.04 


93 


0.13559 


238.88 


94 


0.15735 


241.34 


95 


0.16597 


243.04 


96 


0.17589 


245.52 


97 


0.18018 


246.84 


98 


0.19568 


250.18 


99 


0.16313 


254.00 


100 


0.17123 


257.80 


101 


0.17123 


261.30 


102 


0.17123 


264.80 


103 


0.17123 


268.30 


104 


0.17123 


271.80 


105 


0.17123 


275.30 


106 


0.17123 


278.80 



Prokofiev's Approximation of af 

We choose not to use in the present work experimen- 
tal fission cross sections directly as they are published 
in the literature. Fig. 1 (kindly provided by Dr. 
Prokofiev) explains well the reason: The point is that 
for intermediate- and high-energy reactions, where our 
codes are supposed to be used, the experimental data 
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on proton-induced fission cross sections are sparse and 
not as precise as for low-energy reactions measured for 
reactor applications. Intermediate- and high-energy ex- 
perimental fission cross sections induced by neutrons, 
pions, and other projectiles are even more sparse than 
the ones measured with protons. As one can see from 
Fig. 1, fission cross sections measured at such energies 
in different experiments differ so significantly from each 
other that it is difficult to use such data in development 
and validation of models and codes, without a special 
analysis of all details of every measurement. Fortu- 
nately, this has been done by Prokofiev [13] so we use 
here his results. Prokofiev spent many years on compil- 
ing proton-induced measured fission cross sections and 
on analyzing the details of each experiment. As a re- 
sult, he divided all measurements into three categories: 
1) the highest, where obtained data are very reliable 
and can be used without any mistrust; 2) high-quality 
data, reliable, but requiring some normalization; 3) 
data of low reliability, that would be better not used. 
Then, using only measurements from the first group 
and data from the second group after a corresponding 
re-normalization, Prokofiev developed systematics for 
proton-induced fission cross sections for all preactinidc 
and actinide nuclei for which he was able to find enough 
data [13, 14]. At our energies, we consider Prokofiev's 
systematics as the most reliable "experimental" fission 
cross sections and prefer to use them to develop and 
test our codes instead of using experimental values pub- 
lished in original publications by different authors. 

For subactinide nuclei from 165 Ho to 209 Bi and in- 
cident proton energies above 70 MeV, Prokofiev pro- 
posed [13] the following universal parameterization for 
the proton-induced fission cross section, <jf(E p ) [mb]: 



<j f (E p ) = Pi{l - exp[-P 3 (P p 
x (l-P 4 lnP p ), 



P2)}} 



(4) 



where E p is the incident proton energy [MeV] and Pi, 
P 2 , P3, and P4 are fitting parameters. P4 was fitted as 

P(7 2 /A\-$° if Z 2 /A< 32.32 

n(Z 1 ' ~ 1 Q 4 ,i + Q^Z 2 /A if Z 2 /A > 32.32 



where fitting parameters Qij are given 
rameters Pi, P2, and P3 were fitted as 



(5) 

in Table 2. Pa- 



\{Z 2 /A) ^cxp[Q hl +Q h2 (Z 2 /A)+Q h3 (Z 2 /A) 2 }. (6) 



Table 2. Parameters Qij in the Pi(Z 2 /A) systematics 
for target nuclei from Ho to Bi [13] 



i 


a 1 
3 — 1 


J - 2 


J - 3 


1 


119.0 


-7.852 


0.1332 


2 


9.979 


-0.1847 





3 


-27.40 


0.6792 





4 


-1.140 


0.0352 






For actinide nuclides from 232 Th to 239 Pu and inci- 
dent proton energies above 20 MeV, Prokofiev found 
[13] P 2 = 12.1, P 3 = 0.111, P 4 = 0.067, and 

P 1 {Z 2 IA) = Pn{l - cxp[-R 13 (Z 2 /A - R 12 )]}, (7) 

where R n = 2572, R 12 = 34.99, and R 13 = 2.069. Nu- 
merical values of all Pj parameters of the nuclear targets 
fitted by Prokofiev together with the energy interval of 
fitting are published in Tab. 4 of Ref. [13]. 

In Rcf. [14], Prokofiev extended his systematics to 
describe fission cross sections of preactinide nuclei from 
197 Au to 209 Bi in the energy region from 35 to 70 MeV 
and to predict fission cross sections for nuclei between 
209 Bi and 232 Th, where not a single data point is avail- 
able at present. It was found [14] that one can approx- 
imate fission cross sections of preactinides between Au 
and Bi at proton energies between 35 and 70 MeV with 
the formula 



a(E p ) = a exp 



(E p - E ) z 



2w 2 



(8) 



where Eq = 76.3 MeV. Parameters w and ctq depend on 
the fissioning system and characterize, respectively, the 
steepness and the absolute scale of the fission excitation 
functions and are approximated as following: 

w(A, Z)=a + b(Z 2 /A) + c 5W gs {A, Z), (9) 

where 6W gs is the shell correction to the ground-state 
mass of the fissioning nucleus calculated using the sys- 
tematics of Myers and Swiatecki [15], and a = —33.667, 
b = 1.5699, and c = 0.30069. Parameter ao was fitted 
as 

-(E b -E ) 2 - 



o- = cr b exp^ — J, (10) 

where E b = 70 MeV and a b = a(E b ) is calculated ac- 
cording to the high-energy systematics given by Eq. (4) . 

To predict fission cross sections for nuclei between 
209 Bi and 232 Th at proton energies above 70 MeV were 
there are no data, it was suggested [14] that parameters 
Pi of Eq. (4) can be found by interpolation of the sys- 
tematics [13] predictions. The logarithmic interpolation 
scheme was chosen [14]: 



In Pj = Cn + C l2 x, 



(11) 



where the constants CV, (i = 1 . . A,j = 1, 2) are calcu- 
lated as following: 

x T h In Pi (x B i )- xbi In Pi(xTh) , 10 x 
Oil — , (T^J 

XTh - XBi 
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Figure 1: Experimental proton-induced fission cross sections of 238 U and nat \J nuclei compiled by Prokofiev (symbols) 
compared with results of his sytematics [13] for these cross sections (line). We thank Dr. Prokofiev for sending us 
this figure. 



\nPi(x T h) ~ \nPi(x B i) , . 

Oi2 = , (.13) 

XTh ~ X B i 

where Pi{x) are predictions of the systematics [13] de- 
scribed by Eqs. (5-7), and indexes "Bi" and "Th" 
denote the 209 Bi+p and 232 Th+p fissioning systems, 
correspondingly. The resulting CV, values are: Cn = 
-27.74, C12 = 0.9906, C 2 i = 25.83, C 22 = -0.65 67, 
C 31 = -45.80, C 32 = 1-227, C 4 i = -10.95, and 
C42 = 0.2320. Bellow, we use values provided by Eqs. 
(4-13) to adjust the calculation of fission cross sections 



in our CEM2k+GEM2 and LAQGSM+GEM2 codes. 

Results 

The main parameters that determine the fission cross 
sections calculated by GEM2 are the level density pa- 
rameter in the fission channel, o/ (or more exactly, the 
ratio af/a n as calculated by Eq. (2)) for preactinides, 
and parameter C(Z) in Eq. (3) for actinides. The sensi- 
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tivity of results to these parameters is much higher than 
to fission barriers used in calculation or other parame- 
ters of the model. Therefore we choose to adjust only 
these two parameters in our merged CEM2k+GEM2 
and LAQGSM+GEM2 codes. We do not change the 
form of systematics (2) and (3) derived by Atchison. 
We only introduce here additional coefficients both to 
df and C(Z), replacing a/ — > C a x a/ in Eq. (2) and 
C(Zi) -> C c x C(Zi) in Eq. (3) and fit C a and C c both 
for CEM2k+GEM2 and LAQGSM+GEM2 codes for all 
nuclei and incident proton energies where Prokofiev's 
systematics apply. No other parameters in GEM2 or 
our CEM2k and LAQGSM were changed. For preac- 
tinides, we had to fit only C a - The values of C a found 
by fitting our results to Prokofiev's predictions are close 
to one and change smoothly with changing the proton 
energy and the charge or mass number of the target. 
Such finding gives us confidence in our procedure, and 
allows us to interpolate or extrapolate the values of C a 
for nuclei and incident proton energies not covered by 
Prokofiev's systematics. For actinides, as described in 
[3, 4], we have to fit both C a and C c . The values of C a 
we find are also very close to one, while the values of 
C c are more varied, but both of them change smoothly 
with the proton energy and Z or A of the target, that 
again allows us to interpolate and extrapolate them for 
nuclei and energies outside Prokofiev's systematics. 

We fixed the fitted values of C a and C c in data 
blocks in our codes and complemented them with rou- 
tines for their interpolation/extrapolation outside the 
region covered by Prokofiev's systematics. We be- 
lieve that such a procedure provides quite a reliable 
fission cross section calculation by our codes, at least 
for proton energies and target-nuclei not too far from 
the ones covered by Prokofiev's systematics. Our re- 
sults by CEM2k+GEM2 for preactinides are shown 
in Fig. 2, and for actinides, in Fig. 3. Results by 
LAQGSM+GEM2 are very similar, almost coinciding 
with the ones shown in Figs. 2 and 3, therefore we do 
not duplicate them here. One can see that after fit- 
ting C a and C c , the fission cross sections calculated by 
our codes reproduce very well all the experimental data 
covered by Prokofiev's systematics. 

To see how this approach works for reactions induced 
by other projectiles, we tested our codes on several reac- 
tions induced by neutrons, pions, and photons, without 
any more changes or fitting. Fig. 4 shows several ex- 
amples of such results. We see that our codes describe 
them from quite well to very well, although experimen- 
tal data on pion-induced fission cross sections are not 
so rich and precise, and it is difficult to draw conclu- 
sions from a comparison to this data. The fact that we 
give such fits to fission induced by other probes gives us 
confidence in the value of the fitting procedure we per- 



formed in our CEM2k+GEM2 and LAQGSM+GEM2 
codes. 
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Figure 2: Comparison of Prokofiev's [13, 14] sytematics of experimental (p,f) cross sections of 165 Ho, 173 Yb, 181 Ta, 

183 W; 186 Rej 19Bp t) 197 AU; 202 H& 205^ 204p b) 206p b; 207p b; 208p bj and 209 Bi nudei (lines ) with our present 

CEM2k+GEM2 calculations (circles). 
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Figure 3: Comparison of Prokofiev's [13] sytematics of experimental (p,f) cross sections of 232 Th, 233 U, 235 U 238 U, 
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Figure 4: Comparison of calculated by the modified here CEM2k+GEM2 code fission cross sections induced by 
neutrons on 197 Au and 238 U, 7r~ on 209 Bi and 238 U, and 7 on 208 Pb and 232 Th with experimental data and results 
by previous versions of CEM (see details and references in [1]), as indicated. Experimental data are from: 1) 11: 
Staples [16], Prokofiev [17]; Shcherbakov [18], Eismont [19]; 2) n~: [20]; 3) 7: MAR91 [21], MAR89 [22], TER92 
[25], TER96 [23], TER98 [24], CET02 [26], CAL80 [27], KAP69 [29], LEP87 [28], VEY73 [30], ZHA86 [31]. 
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